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Abstract 

We show that liquid organic scintillator detectors (e. g., KamLAND and Borexino) 
can measure the 8 B solar neutrino flux by means of the u c charged current inter- 
action with the 13 C nuclei naturally contained in the scintillators. The neutrino 
events can be identified by exploiting the time and space coincidence with the sub- 
sequent decay of the produced 13 N nuclei. We perform a detailed analysis of the 
background in KamLAND, Borexino and in a possible liquid scintillator detector 
at SNOLab, showing that the 8 B solar neutrino signal can be extracted with a rea- 
sonable uncertainty in a few years of data taking. KamLAND should be able to 
extract about 18 solar neutrino events from the already collected data. Prospects 
for gigantic scintillator detectors (such as LENA) are also studied. 
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1 Introduction 



Observations of solar neutrinos [1,2,3,4,5,6,7] have offered the first experi- 
mental evidence in favor of non-standard effects, in particular neutrino flavor 
transitions induced by non-zero neutrino masses and mixings. Solar neutri- 
nos have been detected by radiochemical experiments (i. e., Homestake [1], 
Gallex/GNO [2,3] and SAGE [4]) which give an energy-integrated informa- 
tion on the solar neutrino fluxes, and by real time water Cherenkov detectors 
(i. e., Kamiokande, Super-Kamiokande [5,6] and SNO [7]) which allow to ob- 
serve the spectral distribution of solar neutrino events. However, the detection 
threshold in Cherenkov detectors is limited to about 5 MeV by the radiopurity 
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of the target mass and, as a consequence, only high energy 8 B solar neutrinos 
spectrum has been measured. 

In the next future, liquid organic scintillator detectors, such as KamLAND [8] 
and Borexino [9,10], will be operating with the goal of measuring the low en- 
ergy solar neutrino fluxes, in particular 7 Be, CNO and pep solar neutrinos. 
The KamLAND experiment is a 1 kT detector located in the Kamioka mine 
(Japan), at a depth of 2700 m.w.e. of rock, operating continuously since Jan- 
uary 2002, with the main goal of measuring the flux of the P^s coming from 
all the Japanese nuclear power plants. This experiment has spectacularly con- 
firmed the so-called Large Mixing Angle (LMA) solution to the solar neutrino 
problem (see, e. g., [11] for a recent reanalysis). Borexino is a 0.3 kT liquid 
scintillator detector which is being commissioned at Gran Sasso (Italy), under 
3800 m.w.e. of rock, whose main goal is the measurement of the 7 Be solar neu- 
trino flux. Moreover, it has been recently proposed to realize a ~ 1 kT liquid 
scintillator detector, denominated SNO+, at the SNO site (SNOLab, Canada) 
under 6000 m.w.e. of rock, after the completion of the SNO detector physics 
program [12]. It is also under discussion the possibility to realize a gigantic 
(> 30kT) liquid scintillator detector, the Low Energy Neutrino Astrophysics 
(LENA) detector [13], in the Pyhasalmi mine (Finland) at a depth of 1450 m 
(~ 4000 m.w.e.), although other sites (e. g., underwater in the site of Pylos in 
Greece) have also been proposed. The observation of solar neutrinos in these 
detectors, through v— e elastic scattering, is not a simple task, since neutrino 
events cannot be separated from the background, and it can be accomplished 
only if the detectors contamination will be kept very low [9]. Moreover, only 
mono-energetic sources such as 7 Be or pep neutrinos can be detected, taking 
advantage of the Compton-like shoulder edge produced in the event spectrum. 

In this Letter, we show that organic liquid scintillator detectors can also mea- 
sure the 8 B solar neutrino flux by means of the z/ e charged current interaction 
with the 13 C nuclei naturally contained in the scintillators. The possibility to 
use 13 C as a target for 8 B neutrinos was pointed out in the past by [14,15,16]. 
Here, we propose a technique to tag the solar neutrino events. Namely, we 
propose to identify the signal by looking at the time and space coincidence 
with the decay of the produced 13 N nuclei. We perform a detailed calculation 
of the solar neutrino signal and of the background in KamLAND, Borexino 
and SNO+, showing that these detectors will be able to extract the signal with 
a reasonable uncertainty in a few years of data taking. It should be stressed 
that the proposed technique does not involve any modification of the experi- 
mental setup, since one expects a background-to-signal ratio of the order of 1 
or less even assuming the natural isotopic abundance of 13 C (~ 1%) and the 
contamination levels already reached in the KamLAND detector [17]. 

The Letter is organized as follows. In the next section we discuss the neutrino 
interactions with 13 C. In Sec. Ill we calculate the solar neutrino event rates. In 
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Sec. IV we analyze the space and time coincidence with the decay of the pro- 
duced 13 N nuclei. Background issues are described in Sec. V. Sec. VI presents 
the expected sensitivity for KamLAND, Borexino and SNO+ and prospects 
for LENA. In Sec. VII we draw our conclusions. 



2 Neutrino interactions on 13 C 



The 13 C is a stable isotope of carbon with a natural isotopic abundance / = 
1.07%. A small amount of 13 C is, thus, naturally present in organic liquid 
scintillators and can be used as a target for neutrino detection. The relevant 
detection process in our discussion is the charged current (CC) transition to 
13 N ground state: 

v e + 13 C -> 13 N(gnd) + e" . (1) 



The reaction threshold is Q = 2.22 MeV and, thus, only 8 B solar neutrinos 
are detectable (with a neglible contribution from hep neutrinos). In liquid 
scintillators one observes the electron produced in the final state with a visible 
energy which, neglecting the detector energy resolution, is simply equal to the 
electron kinetic energy T e . The cross section of reaction (1) is known with great 
accuracy, since it can be deduced from the decay time of 13 N. One has [15]: 



v(E„) = ^£ Pc E c F(Z,E c ) 
ml ■ ft 

= 0.2167 x 10~ 43 cm 2 ^^F(Z, E e ) , (2) 

where E e = E v — Q + m c is the electron energy, 1 p c is the electron momentum, 
F(Z = 7, E e ) is the Fermi factor and the ft— value of 13 N decay is experimen- 
tally determined as log(/t/s) exp = 3.667 ± 0.001 [18]. By averaging the cross 
section over the 8 B neutrino spectrum, one obtains (a) = 8.57 x 10~ 43 cm 2 , 
which is about one order of magnitude larger than the cross section of u c e — ► 
z/ e e scattering. 

The peculiarity of process (1) is that it can be monitored by looking for the 
delayed coincidence with the positron emitted in the 13 N decay: 

13 N _^ 13 C + ^ + e + ? (3) 

1 We neglected the small recoil energy of the 13 N nucleus (of the order of few keV). 
In this assumption, one has simply T e = E e — m e = E v — Q . 
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which occurs with ~ 99.8% branching ratio (0.2% of 13 N nuclei undergo elec- 
tron capture) and a decay time r = 862.6 s. In this case, the visible energy 
is the sum of the positron kinetic energy and the energy released in e + e _ an- 
nihilation, so that the delayed events have a continuous energy spectrum in 
the range [1.02,2.22] MeV. Moreover, in the absence of macroscopic motions 
in the detector, the 13 N nucleus essentially does not move from its original 
position. The expected displacement due to recoil and diffusion 2 during the 
decay time r is, indeed, smaller than the typical detector spatial resolution, 
a ~ 10 cm. This means that the prompt event produced by the reaction (1) 
and the delayed event produced by the decay (3) have to be observed essen- 
tially in the same position. This condition, as we will see in the following 
sections, is extremely effective in reducing the background. 

Other interaction channels of low energy neutrinos with 13 C can, in principle, 
be considered. First, we discuss the CC transition to 13 N excited states. For 
solar neutrinos, only the lowest excited state (at 3.51 MeV) could be of practi- 
cal importance. The cross section for this process is about 30% of that for the 
ground state and it is calculated theoretically with an uncertainty at the level 
of 30 — 40% [15]. However, the 13 N*(3.51 MeV) decays almost immediately to 
12 C+p with almost 100% branching ratio [19]. As a consequence, it cannot 
be discriminated by the coincidence with the delayed events (3). The other 
relevant process is the neutral current (NC) transition: 

v x + 13 C -> v x + 13 C* . (4) 



Here, only the excited state 13 C*(3.68MeV) is relevant and the excitation to 
other levels has negligible cross section. The cross section for this process, 
averaged over the 8 B neutrino spectrum, is (onc) = 115 x 10 _43 cm 2 [14,16] 
and it is affected by about 30 — 40% uncertainty. The process is in principle 
very interesting since it can give a measure of the total 8 B flux and can be also 
tagged by a monochromatic 7-ray emission back to the ground state. However, 
the great uncertainty and the low cross section make this process hard to be 
competitive with the NC measurement made by SNO. 

In this Letter, we will mainly focus on the information which can be obtained 
from reaction (1), which, at present, seems more interesting in view of the 
larger and much better known cross section, and of the delayed detection 
tagging with the 13 N decay. 



2 Typical value of diffusion coefficients in liquids are D ~ 10 5 cm 2 s 1 , which 
correspond to an average displacement in the time r = 862.6 s equal to l^s = 
V2Dt ~ 0.1 cm. 
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3 The 8 B solar neutrino signal 



A neutrino of energy E u interacting with 13 C through reaction (1) produces 
an electron with kinetic energy T e = E v — Q (neglecting the 13 N recoil). The 
rate R v of prompt events (per unit mass) produced by 8 B solar neutrinos in 
the energy window [T ei i,T e)2 ] is thus simply given by: 

T e , 2 +Q 

Ru = ni3 C $ J &E V a(E v )\(E v )P w {E v ) , (5) 

T e ,i+Q 



where <£> = 5.79 x 10 6 cm~ 2 s _1 is the boron neutrino flux [20], X(E U ) 3 is 
the boron neutrino spectrum [21], a(E u ) is the interaction cross section and 
P ee (E u ) is the electron neutrino survival probability. In the previous relation, 
ni3 C is the number of 13 C atoms per unit mass, which depends on the scintil- 
lator chemical composition: 

nn c = -J2fk— > ( 6 ) 



where I = 1.07 x 10~ 2 is the isotopic abundance of 13 C, u = 1.661 x 10~ 33 kT 
is the atomic mass unit, fk is the mass fraction of the k-th component into 
the scintillator, Xk the stoichiometric coefficient of carbon in the molecule, 
and /j, k is the molecular mass of the k-th molecule. KamLAND scintillator is 
composed by 80% of dodecane (C12H26, with a molecular mass ji = 170.33) 
and 20% of pseudocumene (CgH 12 , with a molecular mass fi = 120.19), which 
correspond to ni3 C = 4.60 x 10 29 kT -1 . Borexino is, instead, composed by 100% 
of pseudocumene, corresponding to nn c = 4.82 x 10 29 kT -1 . The SNO+ liquid 
scintillator composition has still to be decided. For simplicity, we will assume, 
here and in the following, that it will be the same as in Borexino (i. e., 100% 
pseudocumene). 

In Fig. 1 we show the function q(E u ) oc P ee (E v )\(E v )a(E u ) (normalized to 
unity) which gives the relative contribution of neutrinos of different energies 
to the total signal from reaction (1). The solid line is obtained in the assump- 
tion of an undistorted 8 B neutrino spectrum (which can be intended as the 
non oscillatory scenario or a constant suppression of u c ). The dashed line is 



3 A useful approximation for the 8 B spectrum is the following: 
with x = E/Eq, a = 2.92, = 3.49, and E = 14.8MeV. 
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distorted spectrum (<5m 2 = 7.92 X 1 0~ 5 eV 2 , sin 2 2tf=0.86) 
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Fig. 1. Relative contribution of neutrinos of different energies to the total signal 
from reaction (1). Straight line: undistorted 8 B neutrino spectrum (constant Pee)- 
Dashed line: bra 2 = 7.92 x 1(T 5 eV 2 , sin 29 = 0.86. See the text for more details. 

obtained in the assumption of v e — > flavor transitions for the following os- 
cillation parameters which are the current best fit (in 2v) for the whole solar 
and KamLAND data [11]: 



Sm 2 = 7.92 x 10" 5 eV 2 , 
sin 2 29 = 0.86. (7) 

The electron neutrino survival probability has been calculated taking into 
account the Mikheyev-Smirnov-Wolfenstein (MSW) effect in the Sun (for sim- 



Table 1 

Prompt neutrino event rates, delayed energy window efficiency and observed signal 
event rates for KamLAND and Borexino liquid scintillators. See the text for details. 



Prompt event rate 1 Delayed energy windows 2 Signal event rate 1 ' 3 

[2.8, 16] MeV [2.8,5.5] MeV [1.02, 2.22] MeV [1.3, 2.22] MeV [2.8, 16] MeV [2.8, 5.5] MeV 

KamLAND 23.6 6.3 0.77 12.4 3.3 

Borexino 24.7 6.6 1.0 16.8 4.5 
Counts-kTy -1 . 

2 Fraction of 13 N decay events in the delayed energy window. 
3 Applied cuts: TZ = 3, T = 2. 



plicity we have not considered the oscillations in the Earth matter). 4 Neglect- 
ing the detector energy resolution, the function f(T e ) = g(T c + Q) also gives 
the spectral distribution of solar neutrino events, since detection reaction kine- 
matics implies a one-to-one relation (i. e., T e = E v — Q) between the electron 
and neutrino energies. The event spectrum is, in principle, extremely sensitive 
to a possible deformation of parent solar neutrino spectrum. In particular, the 
differences between the two curves in Fig. 1 directly reflect the behavior of 
the electron neutrino survival probability in LMA scenarios. Namely, the rise 
of the LMA spectrum (dashed line) with respect to the standard case (solid 
line) below E u ~ 7 MeV (T c ~ 5 MeV in terms of electron energy) is due to 
the transition from vacuum averaged neutrino oscillations at small energies 
to purely adiabatic transitions at large energies. The observation of this fea- 
ture would be very important as a final confirmation of matter effect in solar 
neutrino oscillations. However, it will be extremely hard to observe it in the 
present detectors, due to the smallness of the expected event rates. 

In first two columns of Tab. 1, we show the neutrino event rates (given in 
counts-kTy -1 ) expected in KamLAND and Borexino scintillators in the en- 
ergy windows [T ejl , T ej2 ] = [2.8,16] MeV and [2.8,5.5] MeV, assuming the 
oscillation parameters in Eq. (7). The lower bound (2.8 MeV) has been cho- 
sen to reduce the background from U — Th contamination. The upper bound 
5.5 MeV has been chosen to focus on the low energy part of the spectrum, 
which, as explained above, is particularly interesting and, moreover, is practi- 
cally unexplored by Super-Kamiokande and SNO. As we can see, the counting 
rates are of the order of 10 — 20 counts-kTy -1 . Moreover, they will be further 
reduced by the cuts, essential to reduce the background. However, as we will 
see in the next sections, the background levels are extremely low so that it will 
be possible to extract the solar neutrino signal with a reasonable uncertainty. 
We remark that the proposed measure does not require any modification of 
the present experimental set-up and, although with a larger uncertainty, can 
be complementary to those coming from SNO and Super-Kamiokande. 



4 A simple and accurate approximation to calculate the adiabatic MSW survival 
probability can be found in [22]. 
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4 Tagging the events 



In order to reduce the background, we can take advantage of the time and 
space coincidence of neutrino events with the positron emitted in the 13 N 
decay. A candidate prompt event will be tagged as signal only if followed by 
a delayed event in the energy window [1.02, 2.22] MeV, within a time interval 
At = Tt (where r = 862.6 s is the 13 N decay time), and inside a sphere of 
radius r = IZa from the prompt event detection point (where a ~ 10 cm is 
the typical detector spatial resolution, see, e. g., [10]). The signal event rate 
is thus given by: 

S = R u • e(T,TZ) , (8) 



where the global efficiency of the coincidence, e(T,TZ), is determined by the 
combined efficiency for the cut in space, (,(7Z), and in time, r](T): 

e(T,K)=Z(K)-Ti(T) . (9) 



The function rj(T) is simply equal to the probability that the 13 N nucleus 
decays within the time At = Tt: 

r](T) = 1 - exp(-T) . (10) 



The function £(Jl) is instead the probability that the prompt event and the 
delayed event, which are assumed to occur in the same position, are detected 
at a distance smaller than r = TZa. Modelling the detector spatial resolution 
with a gaussian function, one obtains: 5 

+,„s J?dxx 2 exp(-x 2 /4) c fH\ /T / U 2 \ 

= roc _, Jl) Ju\ = «rf hr - exp -— . (11) 



/ °° dx x 2 exp(-x 2 /4) V 2 / V tt V 4 



We remark, that the above equation is valid in the assumption that the dis- 
placement between the point where 13 N is created and the point where it de- 
cays is small with respect to a. For that to happen, the macroscopic motions 



5 If we model the detector spatial resolution with a gaussian function /(x, xo) oc 
exp[— (x — x ) 2 /2(j 2 ], where x is the true position of the event and x is the ob- 
served position of the event, the probability that the prompt event is observed in 
the position x p and the delayed in the position Xd can be cast as dP(x p ,Xd) = 
/(x p ,xo)/(x d ,x ) d 3 x p d 3 x d oc exp -(r + y) 2 /4cr 2 d 3 rd 3 y, where r = x p - x d 
and y = x p + x d — 2x . Actually, only the distance r = |r| is observable, so that, 
integrating r on a sphere of radius TZa and y on R 3 , one obtains Eq. (11). 
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in the liquid scintillator have to be sufficiently slow. This can be achieved, 
for example, by maintaining a small temperature gradient pointing upward 
everywhere in the detector. KamLAND data show that the measured average 
displacement of the diffusive 222 Rn over its 5.5 d mean life is less than 1 m 
[23]. Therefore, the assumption that 13 N nuclides displacement over their 15 
minutes lifetime can be kept smaller than detector resolution seems justified. 

Finally, we consider the possibility that the delayed energy window is reduced 
with respect to the full energy range ([1.02, 2.22] MeV) of 13 N decay spectrum. 
In this case the signal event rate is given by: 

E 2 

S = R v -e(T,K)-B(E 1 ,E 2 ) = R v -e(T,K)- J d£ d x (E d ) , (12) 

Ei 

where B(Ei, E 2 ) is the fraction of decay events in the adopted energy window 
[E 1 ,E 2 \ [x(£) is the normalized 13 N decay spectrum]. 

In last two columns of Tab. 1, we show the signal event rates (given in 
counts-kTy -1 ) expected in the KamLAND and Borexino liquid scintillators, 
after that the efficiency cuts are applied. Here, for illustrative purposes, we 
consider a time cut at T = 2 and a space cut at 1Z — 3, which correspond 
to a global efficiency e(lZ, T) = 0.68. It is clear, however, that the cuts must 
be optimized according to each detector's capabilities and performances. For 
the KamLAND detector, moreover, we restrict the delayed energy window 
to [Ei, £2] = [1.3,2.22] MeV in order to reduce the background from 210 Bi 
originated by the decay of 210 Pb that can be either produced by build-up 
due to 222 Rn contamination in the liquid scintillator or caused by an intrinsic 
impurity. 6 This reduces further the efficiency by a factor 0.77. For Borexino 
(and SNO+), we consider instead the full range [1.02,2.22] MeV, assuming 
that the 210 Bi background contribution will be further reduced, since this is a 
pre-requisite for the observation of sub-MeV solar neutrinos. 7 

As a final result, the expected signal event rates are at the level of 10 — 
20 counts-kTy -1 . In order to observe such low counting rates, one clearly 
needs detectors with sufficiently low background levels (and, of course, efficient 
background rejection). Present detectors, as we shall see in the next section, 
already satisfy this requirement. 



6 At present 210 Pb is an important contamination in KamLAND which will be 
removed by purification to allow solar neutrinos detection. 

7 If the detector has an intrinsic efficiency e p (e^) for the prompt (delayed) window, 
the global efficiency is further reduced by a factor e p • e d . We assume for simplicity 
e P = £d = L 
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Table 2 

Rate of background events in the prompt and delayed windows and rate of fake 
coincidences. 





Prompt cncrj 


ly window(s) 


Delayed energy window 1 




Fake coincidences 1,2 


KamLAND 


[2. 


8, 16] MeV 


[2.8,5.5] MeV 


[1.3, 2.22] MeV 


[2., 


3, 16] MeV 


[2.i 


!, 5.5] MeV 


U - Th 




1168 


1168 


1533 










Cosmogcnics (no 12 B) 




20779 


15654 


261555 










Elastic scattering 




967 


558 


398 










Total 




22914 


17381 


263486 




29.1 




22.1 


Borexino 


[2. 


8, 16] McV 


[2.8,5.5] MeV 


[1.02, 2.22] MeV 


[2., 


3, 16] MeV 


[2.i 


!, 5.5] MeV 


U - Th 




1168 


1168 


2373 










Cosmogcnics (no 12 B) 




2968 


2236 


54800 










Elastic scattering 




1004 


577 


3460 










Total 




5140 


3981 


60632 




1.69 




1.31 


SNO+ 


[2. 


8, 16] McV 


[2.8,5.5] MeV 


[1.02, 2.22] MeV 


[2., 


3, 16] MeV 


[2.i 


!, 5.5] MeV 


U - Th 




1168 


1168 


2373 










Cosmogenics (no 12 B) 




32 


24 


583 










Elastic scattering 




1004 


577 


3640 










Total 




2203 


1768 


6416 




0.08 




0.06 



1 counts-kTy 1 . 

2 Applied cuts: TZ = 3, T =2. 



5 The background 

There are three main sources of background for the proposed measure. These 
are: 

i) Internal background due to U — Th contamination and to contamination 
from long lived radon daughters out of secular equilibrium with 238 U (in 
particular 210 Pb); 

ii) Cosmogenic background due to muon-induced production of radioactive 
nuclides, such as n C, 10 C , etc.; 

Hi) Elastic v— e scattering by solar neutrinos. 

These background sources are well known, so that it is possible to perform a 
detailed analysis of their relevance. This is clearly important, because it allows 
to make a realistic estimate. We remark, however, that the real background 
level will be measured directly by the experiments with great accuracy, being 
the background event rate much larger than the signal event rate both in the 
prompt and delayed energy window (before space and time cuts are applied). 

In Tab. 2, we show the contribution of each background source to the to- 
tal background rates (per unit mass) in KamLAND, Borexino and SNO+. 
We have assumed that the U — Th contamination in the detectors is at the 
10~ 17 g/g level, which correspond to the 'present contamination in the Kam- 
LAND detector [17]. The internal background due to the elastic scattering 
of solar neutrinos on electrons has been evaluated assuming z/ e — > z/ MT flavor 
oscillations with the oscillation parameters given in Eq. (7). Finally, the cos- 
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Table 3 

Depth, residual muon flux and average muon energy in the three undergound 
locations considered in the Letter (see [25] for details). 





Depth (m.w.e.) 


M (m^d- 1 ) 


<£ M > (GcV) 


Kamioka 


2700 


230.4 


285 


Gran Sasso 


3800 


28.8 


320 


SNOlab 


6000 


0.288 


350 



mogenic contribution has been obtained by rescaling the results of [24], which 
are relative to Borexino, also to KamLAND and SNO+. This can be done 
by considering that the muon induced background in a given experiment 
scales as oc ni2 C ^ fl (E fJi ) a , where <3> M is the muon flux at the experimental 
site, (Efj) is the average muon energy, ni2 C is the number of 12 C nuclei per unit 
mass in the scintillator (the 12 C is the most relevant target for muon induced 
radioactive nuclei production in liquid organic scintillators) and a = 0.73. 8 
From the data in Tab. 3 (see [25] for details) we calculate that the cosmogenic 
background is ~ 7 times larger in KamLAND than in Borexino, while is ~ 94 
times lower in SNO+. 

From Tab. 2 we see that cosmogenic background is the dominant component 
for KamLAND and Borexino both in the prompt and delayed window, while 
it give only a minor contribution to the total background in SNO+. It is in- 
teresting to give a closer look at the various cosmogenic background sources. 
The relevant cosmogenic nuclei and their half-lives are reported in Tab. 4. We 
see that only n C and 10 C nuclei have a long lifetime. The background coming 
from the other cosmogenics can be efficiently rejected simply extending the 
muonic veto to few seconds after the muon passage through the detector. How- 
ever, in order to avoid further losses of efficiency due to the veto dead-time, 
we rejected only the background events coming from the (fast) decay of 12 B. 
In the last four rows of Tab. 4 we show the event rates (given in counts-kTy -1 ) 
expected in Borexino [24] for the four energy windows considered in this work. 
We see that the main contribution to cosmogenic background in the prompt 
energy windows come from 10 C, while the dominant source in delayed win- 
dows is provided by n C nuclides. It was recently shown [25] that 11 C-induced 
background can be greatly reduced by a three fold coincidence with the parent 
muon track and the subsequent neutron capture on protons. However, in order 
to be extremely conservative, we have not considered this possibility. 9 

8 The cosmogenic production cross section scales with the energy as a oc with 
a ~ 0.73 [24]. 

9 An additional background not included above is provided by the interaction of 
cosmogenic protons with 13 C, according to 13 C(p,n) 13 N. These events are potentially 
dangerous because, being followed by the decay of the produced 13 N, they cannot 
be discriminated by the coincidence (a rejection is, anyhow, possible by looking at 
the subsequent neutron capture on protons). We have estimated, by MonteCarlo 
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Table 4 

Main cosmogenic nuclei, their half-lives, and event rates (in counts-kTy™ 1 ) in the 
Borexino detector [24] for the relevant energy windows considered in this Letter. 



/?+ emitters f3 emitters 





n C 


10 C 


9 C 


8 B 


12 B 


8 Li 


6 He 


T1/2 (s) 


1218 


19.3 


0.127 


0.770 


0.0202 


0.840 


0.807 


[2.8, 16] MeV 





2021 


275 


389 


2933 


241 


42 


[2.8,5.5] MeV 





2021 


40 


82 


911 


51 


42 


[1.02,2.22] MeV 


52630 


1637 


2 


5 


138 


7 


519 


[1.3,2.22] MeV 


35316 


1637 


2 


5 


119 


6 


399 



In the last two columns of Tab. 2, we show the background in the prompt 
energy windows after the coincidence criteria are applied. This is obtained by 
considering that the probability to have a background event in the delayed 
energy window when space and time cuts are applied is simply equal to the 
average number of delayed background events during the time interval At = 
Tt and inside the spherical volume V = (4/3)7r(7£(x) 3 , being this number 
much smaller than one. 10 The rate of fake coincidences B is thus given by: 



B = (B p B d p) 



(13) 



where B p and are the prompt and delayed background rates (per unit 
mass), p is the liquid scintillator density (equal to p = 0.78 g/cm 3 for the 
KamLAND scintillator and p = 0.88 g/cm 3 for the Borexino scintillator) and 
we have taken 1Z — 3 and T = 2. We see that, despite of the large number 
of background events (several thousands per kTy) in the prompt and delayed 
energy window, the fake coincidences are rare (tenth per kTy in KamLAND, 
few per kTy in Borexino or almost absent in SNO+), and comparable or lower 
than the expected signal. For this reason we think that a measure of the 8 B 
solar neutrino flux is feasible. 



simulations, that this background component is negligible, being the cosmogenic 

proton interaction rate with 13 C of the order ~ 10 -2 kTy -1 in Borexino. 

10 The probability to have at least one background event inside a certain (spatial 

and/or temporal and/or energy) window should be calculated by means of a Pois- 

sonian distribution. However, since the number of background events is very small, 

the probability is practically equal to the average number of events in the considered 

window. 
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6 Expected sensitivity and future prospects 



Since real and fake coincidences are indistinguishable, the number of signal 
events Ns has to be obtained from the difference between the total number of 
observed events Nt and the number of background events Nb: 



N s = S ■ E = N T - N B = N T - B ■ 8 



(14) 



where S is the signal event rate, B is the background event rate and S is the 
total detector exposure. The uncertainty of the number of signal events AN S 
is obtained propagating the error in Eq. (14): 



AN S = y/N T + (AB ■ Sy 



(15) 



where we assumed that the total number of events is affected by a Poissonian 
uncertainty AN? = \JNt and we combined in quadrature the errors. Dividing 
by Ns we obtain the fractional uncertainty SS of the signal rate: 



5S = 



AS 



'1 + r 
S-S 



+ r 2 5B 2 , 



(16) 



where r = B/S is background-to- signal ratio and SB is the fractional uncer- 
tainty of the background rate. 

As anticipated in the previous section, the background will be directly mea- 
sured by the experiments. More precisely, one measures the prompt and the 
delayed background rates and, then, determines the final background rate B 
through Eq. (13). The uncertainty 5B can, thus, be estimated as: 



SB = J5BI + 5BI 




(17) 



where we considered that the fractional uncertainties of the prompt and de- 
layed background rates are given by 5B p( i = [-B Pj d • £] _1 / 2 ,as prescribed by 
Poissonian statistics. By using Eq. (17), Eq. (16) can be cast as: 



5S = 



VS~TS\ 



[1 + r) + 



S_ S_ 
lB p B d 



(18) 



It is immediately evident that the second term in the square root of the above 
expression is always negligible, being the ratios r 2 S/B p ^ of the order of 10 -3 
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Table 5 

Background-to-signal ratio and expected sensitivity for KamLAND, Borexino and 
SNO+ for lkTy of exposure. See the text for details. 



Background-to-signal ratio 1 Expected sensitivity 1 Expected sensitivity (optimized) 2 

[2.8, 16] MeV [2.8, 5.5] MeV [2.8, 16] MeV [2.8, 5.5] MeV [2.8, 16] MeV [2.8, 5.5] MeV 

KamLAND 3 2.35 6.73 51.9% 152.2% 51.2% 145.4% 

Borexino 4 0.100 0.291 25.6% 53.6% 23.7% 52.2% 

SNO+ 4 0.005 0.014 24.4% 47.5% 20.6% 40.9% 
1 TZ = 3 and T = 2. 

2 Space and time cuts arc optimized to minimize 5S in each detector, see the text. 
3 Dclayed energy window: [1.3,2.22] MeV. 
4 Delayed energy window: [1.02,2.22] MeV. 



or less in the various detectors. This means that the contribution of the back- 
ground uncertainty to the total error budget is always negligible: 




5S~\— — . (19) 



This situation is, in principle, favorable. The fractional uncertainty 5S scales 
as 8~ 1 / 2 , and, thus, it is possible to obtain a good sensitivity if the detector 
exposure is large enough. In the following we assume 8 = 1 kTy. 

In the first two columns of Tab. 5 we show the background-to-signal ratio r for 
the three experiments under study in the two energy windows [2.8, 16] MeV 
and [2.8, 5.5] MeV, assuming a space cut TZ = 3 and a time cut T = 2. 
In the third and fourth columns, we show the corresponding sensitivity 5S, 
calculated according to Eq. (19), assuming a total exposure 8 = 1 kTy. In the 
last two columns, we give the minimal values for 5S obtained by choosing the 
optimal values of TZ and T which minimize the quantity (1 + r)/S in each 
experiment. 11 Clearly, the lower is the background level in the detector, the 
larger is the space-time window which has to be considered. In KamLAND 
the optimal sensitivity is obtained with TZ = 2.82 and T = 1.65 in the energy 
window [2.8, 16] MeV and TZ = 2.64 and T = 1.44 in the energy window 
[2.8, 5.5] MeV. Larger cuts are preferable for Borexino and SNO+ (of the 
order of TZ ~ 5 and T ~ 5). However, the choice of cuts is not crucial these 
(less noisy) detectors and the sensitivity does not change dramatically if we 
take tighter cuts. 

In KamLAND, due to the large cosmogenic contribution, the background-to- 
signal ratio is equal to about 2.5 in the energy window [2.8, 16] MeV (while 
it is equal to about 7 if we restrict to [2.8,5.5] MeV). This corresponds to 
an expected sensitivity 5S equal to about 50% in one year of data taking 
(assuming ~1 kT fiducial mass). We remind that KamLAND has already 

11 This exercise have to be done numerically, since TZ and T enter in Eq. (19) in a 
non-trivial way. 
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analyzed data for a total exposure £ = 0.766 kTy [8], corresponding to ~ 18 
solar neutrino events in the window [2.8, 16] MeV which can be extracted 
with about 60% uncertainty. 12 Despite the large uncertainty, we believe that 
this measure would, anyhow, represent a milestone, since it would be the first 
observation of solar neutrinos into liquid scintillator detectors. We hope that 
the KamLAND collaboration will try to extract this piece of information from 
their own set of data. 

In Borexino and in SNO+, due to the larger depth of the experimental sites, 
the background-to-signal ratio is much less than one. The sensitivity is thus 
only limited by the statistical error of the signal events. The low background 
level allows to explore with sufficient accuracy the energy windows [2.8, 5.5] MeV 
for which, at present, we have no direct information. This also indicates that, 
in these experiments, it will be possible to decrease the lower bound of the 
energy windows (2.8 MeV) with only a moderate decrease of the expected 
sensitivity. 

We remark that, even if the background is negligible, the low expected count- 
ing rates do not allow to observe a possible distortion of solar neutrino event 
spectrum, unless the 13 C abundance is enriched and/or one considers gigan- 
tic detectors. In principle, 13 C enrichment is possible. 13 However, the current 
separation techniques probably do not allow a massive production of this iso- 
tope. For this reason, we do not consider in detail this possibility. One should 
note, however, that even a partial enrichment (e. g., corresponding to a 13 C 
abundance of the order 20 — 30%) could allow to obtain important results, 
like e. g., the high accuracy determination of the solar neutrino spectrum 
down to energies equal to about E u ~ 3 MeV (or the observation of hep solar 
neutrinos) . 

Finally, we briefly discuss the perspectives for gigantic liquid scintillator detec- 
tors planned in the future. In particular a > 30 kT detector, the Low Energy 
Neutrino Astrophysics (LENA), has been proposed [13]. The site proposed for 
the experiment is Pyhasalmi mine in Finland at a depth comparable to that 
of Borexino (~ 4000 m.w.e.). This means that the cosmogenic background 
will be sufficiently low for the proposed measure. The scintillator should be 
composed mainly by PXE (CieHig, with a molecular mass fi = 210.31 and a 
density density p = 0.998 g/cm 3 ), but, of course, the final composition has yet 
to be decided. It is clear that, in such a large detector, it will be very hard to 
keep the internal background low. However, the gain in statistics will probably 
overcompensate this limitation. To give an example, in one year of data taking 
with a fiducial mass equal to 30 kT, one obtains a better sensitivity than in 



12 Of course, with the proposed cuts, only about 9 of these events would be tagged 
as candidate. 

13 The 13 C is mainly used in health diagnostic, since it has a specific NMR signature. 
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SNO+ (with 1 kT fiducial mass), even assuming a background-to-signal ratio 
of the order of ten. For this reason, we believe that the LENA detector has 
the capability to perform a precise measure of the 8 B flux (comparable to that 
provided by Super-Kamiokande and/or SNO) in a few years of data taking. 



7 Summary and conclusions 

In this Letter we have discussed the possibility to detect 8 B solar neutrinos 
by using the v c CC-interaction with 13 C nuclei naturally contained in or- 
ganic liquid scintillators. The proposed detection process has a low threshold 
(Q = 2.22MeV) and large and well-know cross section. Moreover, one can take 
advantage of the subsequent decay of the produced 13 N nuclei to discriminate 
neutrino events from the background. 

We have calculated the expected event rates (of the order of ~ 20 kTy -1 ) for 
KamLAND, Borexino and an hypothetical Borexino-like experiment situated 
at SNOlab (SNO+). Moreover, we have evaluated thoroughly all the possible 
sources of (external and internal) background in the three considered detec- 
tors. We have shown that the background-to-signal ratio is ~ 2 in KamLAND, 
while is much less than 1 in Borexino and SNO+. 

Finally, we have calculated the expected sensitivity for the various experi- 
ments. Assuming an exposure equal to £ — 1 kTy, the solar neutrino signal 
can be extracted with uncertainty of the order of ~ 50% in KamLAND and 
~ 20 — 25% in Borexino and SNO+. The expected sensitivity scales as £ -1 / 2 , 
since background is directly measured by the experiments. Gigantic (such as 
LENA) and/or enriched detectors, having a much larger statistics, will have 
the possibility to perform a very precise measurement of the 8 B neutrino flux. 

It should be stressed that the proposed measure does not require any modifica- 
tion of the standard experimental set-up. The KamLAND detector should be 
able to extract about ~ 18 8 B solar neutrino events from the already collected 
data (corresponding to a total exposure equal to 0.766 kTy). 
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